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Propagation of sound beam diffracted from a circular aperture in far-field region has been studied 
in this paper by the method of angular spectrum representation and stationary phase method. This 
nonparaxial theory is useful when beam angle is not very small and the wavelength A is comparable 
to aperture diameter a unlike the situation in paraxial approximation. Here we have studied two 
cases, one for a Gaussian source and other for a plane piston source. 

PACS numbers: 43.20.+g,43.35.+d 



I. INTRODUCTION 

Pressure field radiation and diffraction from a source 
with a circular aperture and its propagation through elas- 
tic medium is quite well known and a well studied subject 
in acoustics[2|, |5|, [6( . Paraxial approximation has been 
successfully applied to study such phenomena. This 
approximation neglects term (the second derivative 
of pressure with respect to propagation axial co-ordinate) 
in the Helmholtz equation, as the term is of higher order 
in beam angle 9 which is very small (sin 9 « 9). But 
when 9 is not so small, paraxial approximation is no 
longer valid. Some good works has already been done in 
optics in such nonparaxial correction 0,0, 0- But, in 
acoustics, such work is still missing. This paper concen- 
trates on the propagation of the acoustic beam transmit- 
ted from a monochromatic source from a circular aper- 
ture in nonparaxial far-field zone, through a compara- 
tively general approach employing the angular spectrum 
distribution and the stationary phase method. Gaussian 
and plane piston are taken as the initial source. The pa- 
per begins with a discussion of an analytic nonparaxial 
formalism with analytical results. This is followed by a 
presentation of numerical results. Discussions with con- 
clusions are the final feature of this paper. 



II. NONPARAXIAL THEORY 

A. Helmholtz equation and angular spectrum 
representation 

The governing equation for spatial part of pressure field 
is given by the Helmholtz equation 



[V 2 + k 2 ]P(x,y,z) = 



(1) 



In any plane z = constant, let us assume the field can be 
represented by a Fourier integral as 



P(x,y,z) 
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which when substituted in Eq.l yields 

[V 2 + k 2 ]V{u, v; z)e lk(ux+vy) dudv = (3) 
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and then differentiating under integral sign, we find that 

°° f°° L 2 2 , 2 v«, \ d 2 V(u,v;z)~ 
/ (-u 2 -v 2 + k 2 )V(u,v;z) + y \ ' 

-oo J — oo . ®^ 

xe ik{ux+vy) dudv = 0(4) 

For any values of x and y, the term under bracket will 
go to zero. Hence the function V{u,v;z) will satisfy the 
following differential equation 



d 2 T{u,v;z) 
dz 2 



+ w 2 T(u,v;z) = 



(5) 



Where w 2 = 1 — u 2 — v 2 . When u 2 + v 2 < 1, w 



■v(i 



v 2 ). Again, when u 2 + v 2 > 1, then w 



+iyj{u 2 + v 2 - 1). 

The general solution of the partial differential equation 
(5) is 



V(u, v; z) = A(u, v)e lkwz + B(u, v)e 



— ikv 



(6) 



where A(u, v) and B(u, v) are arbitrary functions of u 
and v. On substitution of equation (2), we obtained 



P(x,y,z) 
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A(u, v)e lk{ux+vy+wz) dudv + 
B(u, v)e- ik(ux+vy+wz) dudv (7) 



— oo J — oo 



If we don't consider reflection mode, then B(u,v) = 
(soft baffle case). Otherwise for rigid baffle, the pressure 
will be doubled due to reflection. So 



V{u,v;z)e lk(ux+vy) dudv (2) P(x,y,z) = 



oo J — OO 
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A{u,v)e lki - ux+vy+wz Uudv, (8) 
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here A(u, v) is defined as: 
k 



A(u,v) 



2- 
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P(x,y,0)e-' kiux+vy) dxdy 
(9) 



Let, the source function be defined as 

P(x,y,0) = P (x,y)T(x,y) 



(10) 



where T(x, y) is the window function for circular aperture 
defined as 



= f 1 if x 2 + y 2 < a 2 
^ ' ' [0 otherwise 

Therefore A(u, v) can be determined by 



(11) 



Po{x,y)e~ lk(ux+vv) dxdy (12) 



x 2 +y 2 <a 2 



Transforming from Cartesian coordinate to cylindrical 
coordinate system i.e. x = pcosO and y — psinO, and 
assuming the source is axially symmetric (Fig.l), the in- 
tegral becomes 



A(u,v) 
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using the identity 

/.2tt 



P (p, Q) e -^ P (ucos0+vs m e) pdf}de 

(13) 



-ikRpcos^-i)^ = 27T J (kRp) (14) 







our problem. Now let us specify some other variables 

s x = x/r, s y = y/r, s z = z/r, 



Jo(fepv u 2 + v 2 ) 

u 2 -\-v 2 <a 2 
xe *kr(us x +vs v +w Sz ) dudv 



u 2 -\-v 2 <a 2 



(18) 



a{u,v)e lkr9{u ' v '' s ^ s ^dudv (19) 



where r = \J x 2 + y 2 + z 2 , a(u, v) = Jo(kp\/ u 2 + v 2 ) and 
g(u,v;s ) = s^u + Syt; + s 2 w. The critical stationary 
points for this integral will be u\ = s x , v\ = Sj, and 
Wi = s z . Hence in this approxination one obtains [J 



2m 



kr^\A 
2m (z 



:a(u 1 ,v 1 )e ikr ^ Ul ' Vl ' Sx ' s y^ 



2m ( z\ Jx 2 + y 2 . e un " 

re \r/ r r 



where 



A — (guugw 5ud)«1:«1 



(21) 



■9" = If >S« = If) 5«« = §Jr etc. Therefore equation ( 16) 
becomes 



k f z\ e 
P{x,y,z) = - I - J 



V* 2 + : 



P (p,0)J (fcp- 
r Jo r 



-)pdp, 
(22) 
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A(«,«) = ^/ P (p,0)J (fcpV / u 2 + V 2 )pd /9 (15) 
27T Jo 



Thus P(x, j/, z) is 

fc 2 '° 



p< •'■■//■ --) - yy p ( (0 ,o)Jo(w« 2 +« 2 ) 

•u 2 +d 2 <a 2 
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pdpdudv 



(16) 



III. PRESSURE FIELDS FOR DIFFERENT 
AXISSYMMETRIC SOURCES 

A. Gaussian source 

Sources with Gaussian amplitude is often been used 
because of its functional form and invariant nature in 
different integral transforms. Let, the source function be 



B. Method of stationary phase 



Po(p,0) =Poe 



(23) 



The integral over u and v on the right hand side of where P = peak amplitude of the source. So, Eq.(22) 

becomes 



Eq.(16) 
I = 



J (kpx/u 2 + v 2 )e lk{ux+vy+wz) dudv 

u 2j rV 2 <a 2 

(17) 

may be approximated by by method of stationary phase. 
Here w = +(1 — u 2 — v 2 ) 2 as u 2 + v 2 < 1, k = — . 
As the evanescent or exponentially decaying wave does 
not give any significant contribution in far field zone, we 
have avoided another condition m — i(u 2 + v 2 — 1) in 



Pok 



P{x,y,z) = — - — 



Jo (a;) can be expanded as 



r J (kp — X +V )pdp 
r 

(24) 



Jo(z)=£(-l) 



m V 4 



(rn 



\\2 



(25) 
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FIG. 1: Schematic diagram of the geometry of circularly aper- 

tured beam at Z = in Cartesian co-ordinate system. ' ' 

line for Gaussian source, line for piston source. 



Thus, we can getQ: 



i f z\ { pp^ r 



xJ2(-l) n (ka 

n=0 

[7(1 +n, 1) - n\ 



2n+2 / x2 ~j~ II 2 



4r 2 



(n 



l\2 



(26) 




FIG. 2: Distribution 



P(x,0,WZ R ) |2 



of a Gaussian sound 



1 P raaI (i,0,10Z R ) 1 

beam diffracted from circular aperture of diameter 3.175mm. 
2(a) frequency=5 x 10 5 Hz, 2(b) frequency=5 x 10 6 Hz 



where 7(0:, /3) is incomplete Gamma function defined as 

7 (a,/3) = / t^e^dt 
Jo 

= a^/Te^ 1^(1; 1 + «;/?) 

where i-Fi(l; 1 + a; /3) is the confluent hypergeometric 
function of the first kind. Eq.(26) is the basic analytical 
result for far- field behavior of nonparaxial Gaussian beam 
with cirlcular aperture. Paraxial result can be obtained 
from this result by expanding r and retaining the first 
term in amplitude part, and retaining up-to the second 
term in phase part, 



B. Piston source 



In piston sources, the amplitude distribution can be 
described as 



P(x,y,0) = P T(x,y) 



(29) 



Here the window function is same as Eq.(ll). For this 
case, the resulting amplitude takes the integral form 



z + 



x 2 + y 2 



2z 



(27) P(x,y,z) 



So that, Eq.(26) simplifies to 



'ikz 



'x 2 +y 2 \ m [7(1 + n, !)-«!] 



X 



Az 2 



(nl) 



H2 



(28) 



Eq.(28) is the Fraunhofer diffraction formula for aper- 
tured Gaussian beam with a circular aperture in paraxial 
regime. 



- — 

r ) r 



J (kp 



\Jx 2 + y 1 



)pdp 



akyjx 2 + y 2 



x J 1 



ak\J x 2 + y 2 



(30) 



This is another result derived in this paper for pressure 
field of nonparaxial circular apertured piston source in 
far-field region. Again, if we expand r in series as given in 
Eq. (27) , we can arrive at paraxially approximated result 
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FIG. 3: Distribution p^f^oiozn) 1 2 °f a plane piston sound 
beam diffracted from circular aperture of diameter 3.175mm. 
3(a) frequency=5 x 10 5 Hz, 3(b) frequency=5 x 10 6 Hz 




Matlab 7.2 co 

IV. NUMERICAL RESULTS 

des has been written to solve Eq.(26), Eq.(28), Eq.(30) 
and Eq.(31). We have taken two different source frequen- 



cies, these are /i = 5 x 10 5 Hz and / 2 = 5 x 10 e Hz. 
Aperture diameter has been taken as 3.175mm. Velocity 
of sound c has been taken as 1500 m/sec. The axial dis- 
tances that has been taken to compute the field is IOZr, 

2 

where Zr is Rayleigh distance defined as Zr = The 
axial distances are 105.5641mm and 1055.6mm for f\ 
and ji respectively. Results has been shown in Fig. 2 and 
Fig.3. 

V. DISCUSSION 

An analytical method has been proposed to investigate 
sound beams beyond paraxial regime. Angular spectrum 
representation and method of stationary phase are the 
two mathematical tools used. The analytical study and 
final results are for a general situation in far-field zone. 
Paraxial results can be obtained as the special cases of 
those results. Eq.(26) is the general nonparaxial solu- 
tion for Gaussian source whereas Eq.(28) is its paraxial 
form. Similarly for the case of plane piston type source, 
Eq.(30) is the general nonparaxial solution and Eq.(31) 
is its paraxial form. Fig. (2) and Fig. (3) are the plots 
of the comparative numerical results for Gaussian and 
plane piston sources respectively. In Fig. 2 (a), nonparax- 
ial and paraxial results differ. Although the maxima is 
same for nonparaxial and paraxial values, but the spread 
is different. But as f increases [Fig. 2(b)], this difference 
diminishes. Similar result has been shown in Fig(3), i.e 
the results for plane piston source. We conclude that the 
analytical results presented here represents a generaliza- 
tion of propagation of sound beam diffracted from cir- 
cular aperture in far-field zone. Not only generalization, 
but a significant correction over paraxial approximation 
too. 
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